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X> Abstract 

Ph We explore the potential of the electric quadrupole transitions 7s 2 S 1 / 2 - 6d 2 -D 3 / 2 , 6d 2 -D 5 / 2 in 

^ radium isotopes as single-ion optical frequency standards. The frequency shifts of the clock transi- 

tions due to external fields and the corresponding uncertainties are calculated. Several competitive 
A Ra, + candidates with A = 223 - 229 are identified. In particular, we show that the transition 
7s 2 S 1/2 (F = 2,m F = 0) - 6d 2 D 3/2 (F = 0,m F = 0) at 828 nm in 223 Ra + , with no linear Zeeman 
and electric quadrupole shifts, stands out as a relatively simple case, which could be exploited as 
O a compact, robust, and low-cost atomic clock operating at a fractional frequency uncertainty of 

10~ 17 . With more experimental effort, the 223 > 225 ' 226 R a + c i oc k s CO uld be pushed to a projected 
Qh performance reaching the 10~ 18 level. 
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I. INTRODUCTION 



Optical atomic clocks based on ultranarrow optical transitions in single laser-cooled 
trapped ions have demonstrated a stability and accuracy significantly better than the 133 Cs 
atom microwave frequency standard. Transitions in various ions are presently under inves- 
tigation as candidates for optical frequency standards, including electric quadrupole transi- 
tions in 40 Ca+ [HE], 199 Hg+ (3H5J, 88 Sr+ [SI E], and 171 Yb+ 019], hyperfine-induced electric 
dipole transitions in 27 A1 + [TU14T2] . and 115 In + [13] and an electric octupole transition in 
171 Yb + [T3]; proposals also exist for 137 Ba + [15] and 43 Ca + [16]. These ion clocks currently 
operate at fractional frequency uncertainties ranging from 10~ 16 to below 1(T 17 , with 
projected accuracies reaching the 10 -18 level. The ultimate performance of each clock de- 
pends on the atomic structure of the ion, the sensitivity of the transition to the external 
environment, and the complexity of the experimental setup needed to operate the clock. 

At our institute an experiment is in progress p2] to measure atomic parity violation 
in single Ra + ions [18J. This experimental setup can be adapted for an investigation of a 
single-ion Ra + clock. In this paper we explore the feasibility of using the strongly forbidden 
electric quadrupole transitions 7s 2 Si/ 2 - 6d 2 D 3 / 2 at 828 nm and 7s 2 Si/ 2 - 6d 2 D 5 / 2 at 728 
nm in a single laser-cooled and trapped Ra + ion as a stable and accurate frequency standard 
[T9H2T] . Our studies are based on the available experimental information about the Ra + 
ion and on many-body atomic theory. The relevant energy levels of 223 > 225 > 226 R a + an d the 
proposed clock transitions are shown in Fig. [TJ The Qd 2 D 3 / 2 and Qd 2 D 5 / 2 levels have a 
lifetime of 600 and 300 ms [20], respectively, corresponding to a Q factor ~ 10 15 for the 
clock transitions. 

A major advantage of Ra + is that all the required wavelengths for cooling and repumping 
and for the clock transition can easily be made with off-the-shelf available semiconductor 
diode lasers, which makes the setup compact, robust, and low-cost compared to clocks that 
operate in the ultraviolet. Moreover, in odd radium isotopes clock transitions are available 
that are insensitive to electric quadrupole shifts of the metastable Qd 2 Dj levels. Such shifts 
are an important limiting factor for several other ion clocks [22]. The radium isotopes under 
consideration are mostly readily available from low-activity sources. 

Optical clocks are important tools to test the fundamental theories of physics. They 
are particularly useful in laboratory searches for possible spatial and temporal variations of 
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FIG. 1. (Color online) 223 > 225 > 226 R a + level scheme with wavelengths taken from Ref. [23J and 
lifetimes from Ref. |20j . The clock transitions are indicated; in 225 Ra + and 226 Ra + two clock 
transitions are considered. 

the physical constants that define these theories. Such searches are strongly motivated by 
cosmological theories that unify gravity and particle physics (see e.g. Ref. [21] )• Laboratory 
tests have placed strong limits on the temporal variation of the electron-to-proton mass 
ratio m e /m p [2"oH2"T] and the fine-structure constant a. The most stringent limit on the 
latter was obtained by comparing two ultrasensitive ion clocks ( 27 A1 + and 199 Hg + ) over the 
period of a year, yielding a limit a/a = (— 1.6 ± 2.3) x lCT 17 /y [25]. The sensitivity to 
a/a results from relativistic contributions to the energy levels that are of order 0(Z 2 a 2 ), 
favoring heavy atomic systems like 199 Hg + . The Ra + clock transition has a comparably high 
intrinsic sensitivity [19j|20l[28] but of opposite sign to that of 199 Hg + , making it a promising 
alternative candidate for testing the time variation of a. Ra + is also very sensitive to 
variations in the quark masses [23 ED] ■ 
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TABLE I. Long-lived neutron-rich isotopes of radium with their lifetime and nuclear spin / |31] . 
magnetic moments /x/ in units of /j, n [32J and quadrupole moments Q in barn [33 j . Also shown are 
the decay series the isotopes occur in, and possible low-activity production sources; A = 227 and 
229 have to be produced by nuclear reactions. 



A 


Half-life 


/ 


W 


Q 


a 




Decay series 


Source 


223 


11.43 d 


3/2 


0.2705(19) 


1.254(66) 


235TJ 


227 Ac (21.8 y) 


224 


3.66 d 













232 rnu 


228 Th (1.9 y) 


225 


14.9 d 


1/2 


-0.7338(5) 







233^j 


229 Th (7.34 ky) 


226 


1.6 ky 













238TJ 


226 Ra, 230 Th (75.4 ky) 


227 


42.2 m 


3/2 


-0.4038(24) 


1.58 


(11) 






228 


5.75 y 













232 rp^ 


228 Ra 


229 


4.0 m 


5/2 


0.5025(27) 


3.09(19) 







The uncertainties were obtained by adding in quadrature the uncertainties given in Ref. (33] . 



II. RADIUM ISOTOPES 

Radium offers a wide range of short- and long-lived isotopes with even and odd nuclear 
spin that could be considered for use as optical frequency standards. Only trace quantities 
of radium are needed to operate a single-ion Ra + clock, but demands on the half-life and the 
ease of production limit the options. The half-life of the isotope should be long compared to 
the excited 6d 2 Dj level coherence time (~ seconds) required to address the ion with lasers. 
Further, it is preferable from an experimental point of view to be able to trap the ions for 
a longer time, at least a few minutes. 

The light (neutron-poor) isotopes A = 209-214, with half-lives that range from several 
seconds up to a few minutes, have been produced at the KVI by fusion-evaporation reac- 
tions [T71 US]- A possible clock candidate could be A = 213, which has a half-life of 2.7 m; it 
is similar to the isotope A = 225, which we consider in detail below. We focus in this paper 
on the heavier (neutron-rich) isotopes with A = 223-229, because they have a half-life of 
longer than one minute, and, moreover, most of them occur in the decay series of uranium 
or thorium and therefore can be produced in sufficient quantities with a low-activity source, 
so that no accelerator is required. Table [I] gives an overview of these isotopes, with their 
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half-life, nuclear spin, and possible production methods. The nuclear magnetic moments 
and quadrupole moments listed are used to calculate the hyperfme constants of the 6d 2 D 3 / 2 
and 6d 2 D 5 / 2 levels of the odd isotopes for which no experimental results are available. 

For Ra + optical-clock purposes, the even isotopes A = 224, 226, and 228, with zero 
nuclear spin, are very similar and spectroscopically relatively simple. They are analogous 
to the 40 Ca + and 88 Sr + clocks. 226 Ra and 228 Ra are available as a source; 226 Ra + can also 
be taken from a 230 Th source, in which case there is no need to ionize the atoms. We 
limit ourselves to 226 Ra + , which is the most easily available isotope, and we consider two 
transitions, namely 7s 2 Si/ 2 - Qd 2 D 3 / 2 and 7s 2 S\/ 2 - Gd 2 D 5 / 2 , as indicated in Fig. [lj 

In the odd isotopes, with nonzero nuclear spin, the presence of hyperfme structure gives 
two advantages. First, in all odd isotopes m F = -H- m' F = transitions exist, which are 
insensitive to the linear Zeeman shift. Moreover, the odd isotopes offer several transitions 
between specific hyperfme levels that in first order do not suffer from the Stark shift due to 
the electric quadrupole moment of the 6d 2 Dj level. In particular, we study the transition 
7s 2 Si/ 2 (F = 2 } m F = 0) - 6d 2 D 3 / 2 (F = 0,mp = 0) in 223 Ra + (no linear Zeeman and 
quadrupole shifts) and 7s 2 Si/ 2 (F = l,m F = 0) - Qd 2 D 5/2 (F = 3,m F = ±2) in 225 Ra + (no 
quadrupole shift), see Fig. [TJ In addition, we consider the transition 7s 2 S\/ 2 (F = 0, m F = 0) 

- 6d 2 D 5 / 2 (F = 2,mp = 0) in 225 Ra + (no linear Zeeman shift), which resembles the 199 Hg + 
clock. We also include the isotopes A = 227 and 229, although their half-lives are rather short 
and they must be produced in nuclear reactions. Specifically, we consider the transitions 
7s 2 S 1/2 (F = 2, m F = 0) - 6d 2 D 3/2 (F = 0, m F = 0) in 227 Ra+ and 7s 2 S 1/2 (F = 2,m F = 0) 

- &d 2 D^/ 2 {F = 0,rriF = 0) in 229 Ra + ; both transitions are free from linear Zeeman and 
quadrupole shifts. 

III. SENSITIVITY TO EXTERNAL FIELD SHIFTS 

All proposed optical frequency standards are sensitive to external perturbations due to 
the electric and magnetic fields present in the trap. These perturbations cause unwanted 
systematic shifts of the frequency of the clock transition. Although for a large part these 
shifts themselves can be corrected for, there is a remaining uncertainty associated with each 
shift due to limited experimental or theoretical accuracy. In this Section, we will investigate 
the sensitivity to the external fields of the candidate Ra + clock transitions for the different 
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isotopes. Input for the required atomic-structure quantities is taken from the recent KVI 
experiment [17] and from experiments at the ISOLDE facility at CERN [321 GHH EH]- The 
wavelengths of the relevant transition are taken from Ref. [23] • When no experimental data 
are available, we rely on atomic many-body theory calculations. 

In the following, we briefly discuss the relevant shifts point-wise. The shift of the clock 
transition is defined as the shift of the excited 6d 2 Dj level minus the shift of the 7s 2 Si /2 
ground state. The results of our calculations for the different Ra + isotopes are summarized 



below and divided into a sensitivity, see Table IV, and an uncertainty, see Table [V} The 
theoretical expressions for the various external-field shifts can either be found in the liter- 
ature or they are straightforward to derive; for completeness, the most important ones are 
given. In the following, we assume that one single laser-cooled radium ion is trapped in a 
radiofrequency (RF) electric quadrupole field, i.e. in a Paul trap. 



A. Doppler shifts 



The motion of an ion in a Paul trap can be described by a secular oscillation with a 
superimposed micromotion oscillation [36] . The micromotion oscillation is directly driven 
by the RF field applied to the trap. Any movement of the ion in the trap can, via the Doppler 
effect, cause broadening and shifts of the frequency of the clock transition. This effect is 
important even when the ion is laser-cooled to the Doppler limit. In the Lamb-Dicke regime 
[37] . which can be reached by Doppler cooling on the strong 7s 2 Si/2-7p 2 Pi/ 2 transition 
at 468 nm, the oscillation amplitude is small compared to the laser-light wavelength, and 
first-order Doppler shifts are essentially negligible [3B1 EH] • Second-order Doppler shifts are 
still present. However, it can be shown that for a heavy ion like Ra + this shift is negligible 
in the Doppler cooling limit [IE], with a projected fractional frequency uncertainty in the 
low 10~ 19 levels. It is, of course, a major challenge to achieve this limit experimentally [TT] ; 
excess micromotion of the ion, caused by electric fields that displace the ion from the middle 
of the RF pseudopotential, needs to be minimized. 
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TABLE II. The available experimental and theoretical hyperfine structure constants (in MHz) of 
the ls 2 Si/ 2 i 6d 2 D 3 / 2 , and 6d 2 D 5 / 2 levels of the relevant odd isotopes of Ra + . The values A'j for 
the isotopes for which no data was available were calculated with A'j = (I/I') x (/•*/ / fJ-)Aj, while for 
Bj we used B'j = (Q' /Q)Bj. The reference values are printed bold. For the 7 s 2 Sij2 two different 
sets of experimental data were available; we used the underlined values. The experimental value 
for Ad of the 6d 2 D 3 / 2 level of 213 Ra was used to calculate Ad for the 6d 2 D 3 /2 levels of the heavy 
isotopes; the 213 Ra magnetic moment used is [ii = 0.6133(18) |32j . There are no data for the Bp 
coefficient of the 6d 2 D s / 2 level, nor for Ad and Bd of the 6d 2 D 5 / 2 level. Consequently, we used 
the theoretical values listed and estimated the uncertainty of the Ad coefficients of the 6d 2 -D 5 / 2 
to be 3%, and the uncertainty of all Bd coefficients conservatively as 10%. 









7s 2 S 1/2 
A s 


Qd 2 D 3/2 
A D B d 


Qd 2 D 5/2 
A D B d 


213 Ra+ 


Expt. 


[35] 


22920.0(6.0) 












Expt. 


m 




528(5) 








223 Ra+ 


Expt. 


[33] 


3404.0(1.9) 










Expt. 


[35] 


3398.3(2.9) 










Theory 


m 


3567.26 


77.08 


383.88 


-23.90 477.09 




Theory 


m 


3450 


79.56 




-24.08 


225 Ra+ 


Expt. 


[33] 


-27731(13) 












Expt. 


[35] 


-27684(13) 












Theory 


m 


-28977.76 


-626.13 





194.15 


227 Ra+ 


Expt. 


[35] 


-5063.5(3.1) 








229 Ra+ 


Expt. 


[35] 


3789.7(2.3) 









B. Zeeman shifts 

Magnetic fields in the trap lead to frequency shifts of the clock transition via the linear 
and quadratic Zeeman effect. For the transitions that suffer from the linear Zeeman effect it 
is hard to quantify the theoretical uncertainty, because the achievable accuracies depend on 
experimental details. In these cases, multiple transitions mp ■<->■ m' F can be used to average 
out the linear effect to the desired level of accuracy. The linear Zeeman shift is absent in 
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Trip = •<->■ m' F = Otransitions, in which case the quadratic Zeeman shift Az/ QZ becomes the 
dominant source of uncertainty. For the state I7, /, J; F, nip) it is given by 



where the magnetic field B is taken along the z-axis; 7 labels all quantum numbers that are 
not specified. We consider only couplings to the hyperfine-structure partners, since other 
contributions will be suppressed; therefore, the quadratic Zeeman effect is negligible in the 
even isotopes. The Zeeman shifts can be calculated from the hyperfine structure constants 
As } d and Bp of the 7s 2 Si/ 2 , 6d 2 D 3 / 2 , and 6d 2 D 5 / 2 levels, and the electron and nuclear 
^-factors. Table [II] lists the available experimental and theoretical values of Ag t D and Bp of 
the relevant odd isotopes. 

1. DC Zeeman shift 

DC Zeeman shifts are caused by the static applied magnetic field present in the trap. 
We assume a magnetic field of 1 mG, which is a typical value needed to split the Zeeman 
degeneracies to order ~10 kHz needed for proper state addressing. Passive shielding of an 
ion trap against magnetic fields has achieved < 10 /zG field stability |41j. This experimental 
number is taken as the uncertainty in the magnetic field strength in Table |V} In order to 
calculate the uncertainty in the resulting shifts, the uncertainties in and Bp, in the 
magnetic field (~ 10 /xG), and in the gj values were taken into account. For gj the free- 
electron values were used with a conservative 1% uncertainty. The uncertainties due to gi 
and the parameters associated with the 7s 2 Si/ 2 state are negligible. 

2. AC Zeeman shift 

The RF voltages applied to the trap electrodes require rather large currents to flow. These 
currents give rise to an AC magnetic field in the trap center. In a perfect geometry, when the 
currents to all electrodes are equal, the individual contributions of the electrodes will cancel 
each other and the net magnetic field will be zero. However, this cancellation could be far 
from complete [25] . The oscillating magnetic field averages over the clock interrogation time 




/iAz/ QZ ( 7 ,J, J,F,m F ) = {gjfi B - g lf i N ) 2 B 2 J(J + 1)(2 J + 1) x 




(1) 
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TABLE III. Dipole scalar, oc\, and tensor, a\, polarizabilities, in units of 47reo«o, ano - quadrupole 
moments, O, in units of eog, for the 7s 2 S\/2-, 6d 2 D 3 /2i an d 6d 2 D 5 /2 levels in Ra + . 



Ref. 



7s 2 S 1/2 



6d 2 D 



3/2 



6d 2 D 5/2 



o --, 
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[21J 



104.54(1.5) 
106.22 



83.71(77) 

-50.23(43) 
2.90(2) 



82.38(70) 

-52.60(45) 
4.45(9) 



(which is of the order of the 6d 2 D j-level lifetime), which is long compared to typical RF 
periods (0.1-1 /is). Therefore, the expressions for the DC Zeeman effect can be used, with a 
rms magnetic field. For the 199 Hg + clock this magnetic field is conservatively estimated to 
be of the order ~ mG [25]. We use 1 mG as estimate in Table [Vj because for Ra + the mass 
and other trapping parameters are similar. The resulting AC Zeeman shift proves to be one 
of the largest shifts. Therefore, it is important to work with a rather weak trap potential, 
as the average magnetic field scales with RF power. By varying the trap parameters the 
AC Zeeman shift can be measured. Moreover, averaging schemes that exploit the hyperfine 
structure could significantly reduce the uncertainty in the AC Zeeman shift. In this way it 
should be possible to reduce this uncertainty to the level of 25% of the shift itself; this is 
the uncertainty used in Tabel |Vl 



C. Stark shifts 



Stark shifts result both from static electric fields (causing DC Stark shifts) and from 
dynamic electric fields (causing AC Stark or light shifts). First, quadratic dipole Stark 
shifts are discussed, which are caused by the interaction of the dipole moment of the atom 
with the electric field. Next, we discuss quadrupole Stark shifts, caused by the interaction 
of the quadrupole moment of the atom with the gradient of the electric trap field; we look 
at both linear and quadratic quadrupole Stark shifts. 
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1. DC dipole Stark shift 



The theory of the static quadratic dipole Stark shift was developed by Angel and Sandars 
[42J. For the state (7; J, m j) this shift is given by 

hAv nCJ3S (l, J, mj) = -iaj( 7 , J)E 2 - ^(7, J) + x) ~ (3ff - , (2) 

where -E 1 is the DC electric field strength, and a\ are the scalar and tensor polarizabilities, 



respectively. In Table |III| the available theoretical calculations for these polarizabilities are 
listed for the 7s 2 Si/ 2 , 6d 2 D 3 / 2 , and 6d 2 D 5 / 2 levels in Ra + ; we used the results of Ref. [2T] in 
our calculations. The polarizabilities for the hyperfine levels I7, 1, J; F, mp) are calculated 
using 

ai( 7 , /, J, F) = (-l) J +'+™(2F + 1) J F F k \ ci(i, J) . (3) 

\J J ij 

For an ion laser-cooled to the Lamb-Dicke regime, DC electric fields at the position of the 
ion can be reduced to < 10 V/m in the process of minimizing the micromotion [25J. This is 
the field uncertainty that we assume to estimate the fractional uncertainty in Table [V] in a 
worst case scenario, i.e. E z = E. 

The main source of DC dipole Stark shifts, however, is the presence of black-body (BB) 
radiation due to the nonzero temperature T of the trap and its surroundings. The energy 
shift of a level with dipole scalar polarizability in a BB electric field is given by [13] 



/»Ai*b( 7) J,mj) = ~\ (8.319 V/cm) 2 (j^J 0^(7, J)(l + v) , (4) 

where r] is a small calculable term associated with dynamical corrections; it is of the order of 
a few percent jH] and therefore it can be neglected compared to the overall 10% uncertainty 



given in Table IV, which is mainly due to the theoretical uncertainties in the polarizabilities. 
The BB radiation is assumed to be isotropic, so the tensor polarizability plays no role. Since 
the BB radiation shift results in a relatively large fractional frequency uncertainty at room 
temperature T = 293 K (see Table |VJ), the calculation was also performed for liquid-nitrogen 
temperature, T = 77 K (the 199 Hg + clock operates at 4 K). We assume an uncertainty in 



the temperature of 1 K, as in Ref. |45]. 
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2. AC dipole Stark shift 



The most important cause of AC dipole Stark shifts is the laser locked to either the 728 
nm or 828 nm clock transition, since we assume that the cooling and probing lasers are fully 
extinguished at the time of measurement. When the laser propagates along the z-axis, the 
AC dipole Stark shift of a state (7, J, mj) is given by [19] 

7 a / t \ II ( 1/ \ a 1/ \ m J 1/ \ 3m 2 , — J( J + 1) \ 

frAz/ ACDS (7, J,mj;u L ) = -— I a Q [y L ) + Aa^) — - a 2 [y L ) gW/^Tj J ' ^ 

where II is the intensity of the laser which we take as 1 /iW/mm 2 , Ul is its frequency at the 
clock transition, and A is a numerical factor whose value depends on the type of polarization. 
Further, a^z/^), aj(z/^), and ot\{yL) are the dynamic scalar, vector, and tensor polarizability, 
respectively, of the state (7, J, raj). We choose the polarization such that A = 0, therefore 
we only need the scalar and tensor polarizabilities. These are given by 

9 a p 

alh,J-v L ) = — = KVJ'HDlbJ)! 2 =^ , (6) 

0U ; 3(2J+1)4^ (AE) — {hvi) 2 KJ 

^ J - ^4 ( (2, + 3)( ( f + "l)(2, + l) ) V2 X 

y(-i) j - j 'r 1 2 \\(iA\D\\iJ)\ 2 — ^ — , (?) 

with AE = E — E' and D the dipole operator. For ul — > 0, the above equations reduce 
to their static counterparts. In calculating the dynamic polarizabilities we use the values 
for the dipole matrix elements given in Refs. [101 EH EI]- In using this sum over the 
valence states approach, we do not take the core contributions, which are of order 10% [20J, 
into account. However, the core contributions cancel since we look at differential shifts, and 
these contributions are common. The remaining uncertainty is due to neglected higher-order 
valence and valence-core couplings, and the uncertainty in the dipole matrix elements. 

3. Quadrupole Stark shift 

The interaction of the atomic quadrupole moment with the gradient of an electric field 
gives rise to an electric quadrupole shift. This shift is troublesome in several optical frequency 
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standards [22]. The expression used for the linear quadrupole Stark shift is 

2 [F(F + 1) - 3m%] (2F + 1) 



^Az/ LQS (7, I, J, F, m F ) = v4 DC 9(7, J) 



[(2F + 3)(2F + 2)(2F + 1)2F(2F - 1)] 



1/2 



x(-l) /+J+F <( " H - I X , (8) 





2 






2 




I* 


/ 









I) 



-1 



where A DC is the electric field gradient, 0(7, J) the quadrupole moment, and X contains the 
angular factors resulting from the rotation of the quadrupole field frame to the quantization 
axis [46J. The quadrupole moment of the 7s 2 S'i/2 ground state is zero, those of the 6d 2 D 3 / 2 



and 6d 2 Ds/2 levels (20] are listed in Table III There are three special cases in which the 
first-order effect also vanishes for particular hyperfine states of the 6d 2 Dj levels: 

(i) F = levels have no quadrupole moment; this applies to the 223,227 ' 229 Ra + cases. 



(ii) When F = 2, I = 3/2, J = 3/2, the 6j-symbol in Eq. ([8]) is zero. This set of quantum 
numbers is available in 223,227 Ra + , however, there is no improvement over the previous 
case (i). All other shifts and associated uncertainties were calculated to be equal to, 
or larger than, their counterparts in the F = case. Therefore, these transitions have 
not been included in Tables [IV] and [Vj 

{iii) For F = 3, nip = ±2 the shift vanishes because of the factor F(F + 1) — 3m F in Eq. 
Q; this applies to the 225 Ra + (l) case. 

The transitions in 226 Ra + and 225 Ra + (2) do suffer from a linear quadrupole shift. These 



are given in Tables IV and [V} To estimate the size of these shifts and their uncertainties, 
we assumed that in the trap a typical static stray electric field gradient A DC ~ 10 3 V/cm 2 
is present due to patch potentials. We assume that the angular factor X is of order 1. 
Since the orientation of the stray field is unknown, we take the full shift as an estimate 
of the uncertainty. The effects of the much larger RF trapping fields average out over the 
interrogation period. 

However, the transitions that are free from the linear effect do suffer from a second-order, 
quadratic quadrupole Stark shift. This contribution is significant because now the effects 
from the RF trap potential do not average out. This RF potential gives rise to a typical 
rms field gradient A AC = 10 4 V/cm 2 . To estimate the size of the shift we assume that the 
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magnetic-field orientation and the z-axis of the quadrupole trap field coincide. Taking only 
couplings to hyperfine partners into account results in 



hAu QQS ( 7 , 1, J, F, m F ) = AAl c B(l, J? £ (2F + ^ + 1 




-m_F rnj? 

It should be feasible to achieve an overall 10% accuracy in the determination of this shift, 
which is the uncertainty quoted in Table |V} 

It can be seen in Table |V| that, similar to other clocks, the linear quadrupole shift is by 
far the largest shift in Ra + . In 199 Hg + it was canceled by means of an averaging scheme 
[251 HH1 [52], which brought down the uncertainty level to the 10~ 17 level. An alternative 
was presented more recently for 88 Sr + in Ref. [TJ, where it is projected that the uncertainty 
caused by the electric quadrupole shift can be reduced to the 1CT 18 level. 

IV. DISCUSSION AND CONCLUSIONS 



Tables IV and M contain the quantitative results of our studies. Table IV lists the sensi 



tivities of the isotopes under study to the external fields. Also in Table [TV] the sensitivities 
of three other ion clocks that are based on an electric quadrupole transition are shown for 
comparison. In Table |VJ the Ra + sensitivities have been combined with typical values (and 
uncertainties) for the required and spurious external fields to quantify the resulting shifts 
and the fractional frequency uncertainties Sv/u, where v is the transition frequency and Su 
the uncertainty in the total shift. In the top half of the Table the different shifts are given 
in mHz, with the corresponding uncertainty between brackets. 

The transitions 225 Ra + (l), 226 Ra + (l), and 226 Ra + (2) suffer from the linear Zeeman (LZ) 
shift, which therefore has to be controlled to the desired level of accuracy. The transitions 
225 Ra+(2), 226 Ra+(l), and 226 Ra+(2) suffer from a linear quadrupole Stark (QS) shift of 
the order 6-24 Hz, which has to be cancelled in order for these cases to be competitive. 
As mentioned, an averaging scheme was implemented for 199 Hg + , a system comparable to 
225 Ra + (2), and 1CT 17 levels have been achieved j25j |52]. With an alternative averaging 
scheme, it appears feasible to reduce the QS shift experimentally to the 10~ 18 level in 88 Sr + 
[7J, a system comparable to 226 Ra + . The transitions in 223 Ra + , 227 Ra + , and 229 Ra + are 
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insensitive to both the LZ and the linear QS shifts, which is in principle a clear experimental 
advantage. The quadratic QS shifts are only of the order of 1 mHz. 227 Ra + is overall slightly 
better than 223 Ra + , while 229 Ra + is worse, because it has a relatively large quadratic Zeeman 
(QZ) shift. As discussed, of these three, only 223 Ra can be obtained from a source. 

Provided that the LZ and linear QS shifts can be cancelled in 225 ' 226 Ra + , the largest 
remaining shift is caused by the BB radiation. It is of order 0.2 Hz in all the isotopes. As in 
the case of 199 Hg + , this shift can be rendered negligible by cooling down the system, albeit at 
the cost of a more complicated experimental setup. For that reason the BB shift is given for 
two temperatures, namely for room temperature (293 K) and for liquid-nitrogen temperature 
(77 K). The combination of these two options with the possibility of averaging away the QS 



shift (indicated by "no QS" in Table V) give us in total four different results for four sets 
of experimental choices, as shown in the bottom half of Table |Vj In the calculation of these 
uncertainties in the case "no QS," we have assumed that the LZ shift and the linear QS shift 
can be averaged out experimentally to negligible values. The actual obtainable accuracies 
in these cases depend on experimental details, but, as discussed, it appears realistic to aim 
for accuracy levels of a few times 10~ 18 . 

We conclude that in particular the isotopes 223 > 225 > 226 R a + are promising clock candidates 
with projected sensitivities that are all below the 10~ 17 level. The actual experimental 
feasibility of the scenarios discussed above remains to be demonstrated, of course. 223 Ra + 
stands out as an attractive simple candidate, without LZ and linear QS shifts, providing a 
compact, robust, and low-cost atomic clock. 



V. SUMMARY 



In summary, a theoretical analysis of the possible performance of a radium single-ion 
optical clock was presented. It was shown that transitions in several readily available Ra + 
isotopes are excellent candidates for alternative optical frequency standards. The advantages 
of a heavy single ion that can be directly laser-cooled and interrogated with off-the-shelf 
available semiconductor lasers are clear for many applications in which costs and system 
size and stability are of importance. Furthermore, Ra + is an excellent laboratory for the 
search for variation of fundamental constants, where it ranks among the most sensitive 
candidates. 
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